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By exploiting density-matrix renormalization group techniques, we investigate the dynamical spin structure factor of 
a spin-1/2 Heisenberg chain with ferromagnetic nearest-neighbor and antiferromagnetic next-nearest-neighbor exchange 
interactions in an applied magnetic field. In a field-induced spin nematic regime, we find gapless longitudinal and gapped 
transverse spin excitation spectra, in accordance with quasi-long-ranged longitudinal and short-ranged transverse spin 
correlations, respectively. The gapless point coincides with the dominant longitudinal spin correlation, whereas the gap 
position exhibits a characteristic field dependence contradicting the dominant transverse spin correlation. 


Frustrated quantum magnets are a class of interacting spin 
systems in which there are competing interactions that cannot 
be satisfied simultaneously. In general, the combined effects 
of frustration and quantum fluctuations prevent conventional 
magnetic order, and we envisage the emergence of novel spin 
states. Even when the system is magnetically disordered, we 
typically find “hidden” order described by multiple-spin order 
parameters, such as dimer, chiral, and multipole ones. In this 
paper, we particularly discuss a quadrupole state, called a spin 
nematic (SN) state.A schematic picture of the SN state is 
that spins align in a spontaneously chosen axis, while they still 
fluctuate within the axis and their directions remain unfixed, 
resembling the directional order of rod-shaped molecules in 
the nematic liquid crystal. 

The occurrence of the SN state has been demonstrated for 
several kinds of frustrated quantum spin models.Among 
them, a spin-1/2 chain with ferromagnetic nearest-neighbor 
/i < 0 and antiferromagnetic next-nearest-neighbor /2 > 0 
exchange interactions in a magnetic field h, described by 

// = y, S, • Sm + 72X -hYjS], (1) 

i i i 

has been attracting much current interest because of not only 
the novelty of a theoretically predicted SN state,but also 
the relevance to a wide variety of edge-shared copper-oxide 
chain compounds. 

On the theoretical side, the ground-state phase diagram has 
been studied in detail,^’and it is now well established that a 
SN state exists at high fields, where the quadrupole correlation 
{S f S ^S exhibits a quasi-long-range order due to the 
formation of two-magnon bound states. The longitudinal spin 
correlation (5 ^5^.) also shows a quasi-long-range order, while 
the transverse spin correlation {SfS^} decays exponentially. 
Accordingly, the SN state should have gapless excitations in 
the quadrupole and longitudinal spin channels, whereas there 
should be a finite excitation gap in the transverse spin channel, 
corresponding to the binding energy of a magnon pair. Note 
that in the SN regime, the quadrupole correlation is dominant 
below the saturation, and the longitudinal spin-density-wave 
(SDW) correlation becomes dominant as the magnetic field is 
decreased, although both correlations are quasi-long-ranged. 
Thus, hereafter, we refer to this phase as the SN/SDW phase. 
When the magnetic field is further decreased, the ground-state 
phase changes to a vector chiral (VC) phase at low fields. 


Concerning the experimental realization of the SN state in 
the fruatrated ferromagnetic chain, LiCuV 04 has been studied 
frequently as a prototypical material. In a field-induced 
phase, neutron diffraction experiments have found a cohinear 
spin-modulated structure, which is consistent with theoretical 
results for the longitudinal SDW correlation in the SN/SDW 
phase. The spin correlations exhibit short-range behavior 
in ah directions.These observations are indeed suggestive 
of the development of the quadrupole correlation. However, 
since the direct observation of the non-magnetic SN state is 
difficult, the order parameter is still not identified yet. 

To collect evidence that the SN state occurs in reality, it is 
important to find any specific indications in what we observe 
by using various microscopic probes. In this context, recent 
theoretical studies have pointed out that the NMR relaxation 
rate shows characteristic temperature and field dependencies 
in the SN/SDW phase.NMR experiments for LiCuV 04 
have reported consistent results with theoretical predictions, 
signaling the formation of the bound magnon pairs.The 
properties of the low-energy spin excitation spectra have also 
been discussed.^‘^’^^“^^^ 

In this paper, to clarify the property of the SN/SDW phase 
from the viewpoint of the spin dynamics, we investigate spin 
excitation spectra in a wide range of momentum and energy 
by numerical methods. We clearly find that the longitudinal 
spin excitation spectrum is gapless, while the transverse one 
is gapped, as naively expected from the behavior of the spin 
correlations. We discuss the field dependence of the spectral 
weight transfer in relation with the dominant spin correlation. 
A striking feature is that the momentum position of the gap 
deviates from that of the dominant transverse spin correlation 
as the system approaches the saturation. 

Let us consider the model (1) in an A-site chain, and take 
/2 = 1 as the energy unit. We investigate the spin excitation 
dynamics at zero temperature by exploiting density-matrix 
renormalization group (DMRG) techniques.We employ 
the finite-system algorithm in open boundary conditions. We 
compute the dynamical spin structure factor, defined by 

S^(q,oj) = Aim{^^lSf ^ (2) 

71 ^OJ + Eq-H + VTI^ 

where |j/^g) is the ground state with eigenenergy Eq. Note that 
|j/^g) is given by the lowest-energy state in the subspace of a 
given magnetization m = M/N, where M = For the 
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Eig. 1. (Color online) Intensity plots of the dynamical spin structure factor uj) at /i = -1 and /2 = 1, obtained by dynamical DMRG calculations with 
N = 128. (a) uj) at m = 0 and /z = 0. (b) S^{q, u) and (c) S^{q, oj) at m = 0.125 (= 16/128) and h = 0.649. 


calculation of co) at m, we set the magnetic field to be 
the midpoint of the magnetization plateau of M in the /V-site 
system. 7 / is a small broadening factor, and we set 77 = 0.1 
unless otherwise specified. Our particular interest is to clarify 
the anisotropy between longitudinal S^{q,oj) and transverse 
S^{q, oj) due to the formation of the two-magnon bound states 
in the field-induced SN/SDW phase. 

For the analysis of the spin excitation spectrum, we use a 
dynamical DMRG method,^^^ targeting the ground state \iI/q), 
an excitated state Sq\il/Q), and the so-called correction vector 
[oj Eq - H 'iq\~^Sq\il/Q). Here, we note that the truncation 
error rapidly increases as the number of target states increases. 
Therefore, to obtain S^{q, oj) with keeping high accuracy, we 
calculate S^{q,oj) and S~{q,oj) separately, and then use the 
relation 5 oj) = [S^(q, oj)-\-S ~(q, oj)]/4, instead of directly 
calculating S^(q, oj). Note also that we calculate the spectrum 
at q and oj after one DMRG run with fixed q and oj, so that 
we need to perform a great number of DMRG runs to obtain 
a full spectrum. 

Let us start with a brief discussion on the spectrum at zero 
field. In Fig. 1(a), we present the intensity plot of S^(q, oj) at 
/i = - 1 , /2 = 1 (energy unit throughout the paper), m = 0, 
and h = 0. Note that S^(q, oj) = S^(q, oj) aX h = 0 due to the 
SU(2) spin rotation symmetry. We find a sinusoidal dispersion 
that gives the lower boundary of a continuum. The sinusoidal 
dispersion represents the spinon excitation, described by the 
des Cloizeaux-Pearson mode,^^^ since the system decouples 
into two antiferromagnetic chains if /i ^0. We see a large 
amount of spectral weight at a lowest-energy peak (qo, ojq) = 
(7r/2,0.04). Note here that the spectrum is asymmetric with 
respect to qo, as pointed out by the previous studies.That 
is, the spectral weight mainly lies near the lower boundary of 
the continuum for q < q^, whereas it is distributed to a high 
energy region for q > q^. On the other hand, ojq coincides 
with a spin excitation energy, 

A(/V,M) = [Eq{N,M+\) + Eq{N,M-\)-2Eq{N,M)]I2, (3) 

where Eq{N, M) is the lowest energy of the A^-site system in 
the subspace of M at /z = 0. We find that A(/V, 0) shifts to 
lower energy as N increases, and it is extrapolated to almost 
zero in the limit of N ^ 00 . Note that an exponentially small 
gap has been predicted by renormalization group theory, but 
it is hard to detect such a small gap numerically.^^^ 

Now, let us look into the longitudinal and transverse spin 


excitation spectra in the SN/SDW phase. In Figs. 1(b) and 
1 (c), we show S^{q, oj) nndS^^q, oj), respectively, at /i = - 1 , 
m = 0.125, and h = 0.649. For 5'^(^, oj), a lowest-energy peak 
is at {qo, oj^) = (0.375;:, 0.00). That is, q^ moves toward small 
momentum from the position at zero field, ojq is nearly zero, 
indicating a gapless mode for the longitudinal spin excitation. 
We see that a certain amount of spectral weight is transferred 
to the origin due to a finite uniform magnetization, which is 
also indicative of a gapless mode. These gapless points are 
clearly visible due to the large intensity. Moreover, S^(q,oj) 
seems to be gapless at ^ = ;r - ^0 and q = n, although we 
do not observe significant intensity near the possible gapless 
points in ^ > 7i/2. 

In contrast, for S^(q, oj), we observe a lowest-energy peak 
at (qo, ojq) = (31;r/64,0.15). That is, qo remains near;r/2 even 
in the SN/SDW phase. On the other hand, ojq appears to be a 
finite energy. In fact, ojq agrees with the spin excitation energy 
A(N, M), and it is extrapolated to a finite value in the limit of 
N ^ 00 with m = M/N fixed, indicating a gapped mode for 
the transverse spin excitation. Note that S^(q,oj) consists of 
S~(q, oj) and 5'^(^, oj). As shown in Fig. 2, both spectra have 
the lowest-energy peak at the same position, while the overall 
structure is different between them. We find that S~(q, oj) is 
highly dispersive in a wide range of momentum and energy, 
and S'^(q,oj) is less dispersive and it is mainly concentrated 
in a small region near the lowest-energy peak. 

Here, let us discuss how the lowest-energy peak position 
depends on m in the whole range of m. In Fig. 3(a), we plot 
qo of S^(q, oj) as a function of m for several values of Ji. Note 
that we also find a sharp peak at the origin for finite m, but 
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Eig. 2. (Color online) Intensity plots of the dynamicsl spin structure factors 
(a) S~(q,aj) and (b) S^(q,aj) at Ji = -1, /2 = ^ = 0.125, and h = 0.649. 

Note that S^(q, oj) = oj) + S~(q, uj)]/4 is plotted in Fig. 1(c). 
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Eig. 3. (Color online) The lowest-energy peak position of (a) S^{q,(jj) 
and (b) S~iq, lo) for several values of J\ at /2 = 1 as a function of m. Solid 
symbols denote the VC phase at low fields, and open symbols denote the 
SN/SDW phase at high fields. Note that qo shows a stepwise change simply 
due to finite-size effects. The resolution of the momentum is InjN and we 
use V = 128 in the present calculations. Here, the broadning factor is set to 
q = 0.02 to determine the position of the lowest-energy peak precisely. 


we focus on the field-induced shift of the peak position which 
originally locates near njl at zero field. At small m, where 
the system is in the VC phase, shows little dependence on 
m. At large m, where the system is in the SN/SDW phase, qo 
follows the relation qo = (1/2 - m)7i regardless of /i, which 
supports the bosonization result.^"^^ Note that q^ is represented 
by the density of bound magnons 112- m. 

In Fig. 3(b), we show the m dependence of qo of S~(q, oj). 
In the VC phase at small m, we notice that qo agrees with the 
incommensurability of the transverse spin correlation, as will 
be shown in Fig. 4(d). Note here that it has been revealed that 
the incommensurate wave number in the VC phase is strongly 
quantum renormalized toward 7r/2 compared with the pitch 
angle cos“H-/i/4/2) of the helical order in the classical spin 
case.^^ Indeed, we observe that qo = 7r/2 at Ji = -0.5 and -1, 
although the corresponding classical pitch angle is 0.46;: and 
0.42;:, respectively. For small |/i|, we see that qo remains at 
71/2 below a threshold value of m even in the SN/SDW phase 
at large m. However, qo goes away from 7r/2 as m approaches 
the saturation, and the deviation from 7t/2 is more pronounced 
for larger |/i|. On the other hand, the exact energy dispersion 
of the one-magnon excited state in the fully polarized state is 

€i(q) = Ji(cosq - 1) + J2(cos2q - 1) h, (4) 

and its minimum is at ^ = cos“^(-/i/4/2), which coincides 
with the classical pitch angle.^^ The present numerical results 
at the saturation totally agree with this exact description. We 
mention that the bosonization analysis shows that the bottom 
of the one-magnon band is found at 7r/2,^^^ but it is valid only 
for the weak-coupling regime |/i| J 2 and inapplicable in 
the limit of m ^ 1/2. We thus confirm that the bottom of 
the one-magnon band moves from the quantum renormalized 
incommensurate wave number to the classical pitch angle as 
m increases from zero to the saturation. This feature would be 
useful to determine exchange couplings of real materials. 

To gain an insight into the properties of the field-induced 
spectral weight transfer and the dominant spin correlation, we 
examine the spin structure factor, 

S%q) = 5“|.Ag>, (5) 

with attention to sum rules on the integrated intensity of the 





Eig. 4. (Color online) (a) S^{q) at Ji = -1, (b) S^{q) at /i = -1, and (c) 
S^{q) at J\ = -2.5 for several values of m. (d) Pitch angle Q, determined 
from the maximum point of S~{q), as a function of m. Solid symbols denote 
the VC phase at low fields, and open symbols denote the SN/SDW phase at 
high fields. /2 = 1 is fixed. Data are obtained by DMRG calculations with 
N = 128. 


dynamical spin structure factor. To obtain S^{q), we perform 
ground-state DMRG calculations independent of dynamical 
DMRG runs for 5'^(^, co), since we can obtain accurate results 
with relatively small computational cost. The sum rule for the 
energy-integrated intensity reads 

I\q) = ^^S\q,oj) = S\q), ( 6 ) 

leading to the spin structure factor. As for the total intensity 
‘^^(^)^ we have 

f = F = N/4, I- = N/2 + M. (7) 

In Fig. 4(a), we show S ^(q) for several values of m at /i = -1. 
There is a clear peak at ^ = 7r/2 for m = 0, and it shifts 
toward small momentum with increasing m, as indicated by 
vertical arrows. The peak position of S^(q) representing the 
SDW correlation is in agreement with the lowest-energy peak 
position qo of S^(q,(jj) in Fig. 3(a). We also find that S^(q) 
grows at ^ = 0 as m increases. In other momentum parts, 
S^(q) is suppressed so as to keep the total intensity F = A/4. 
Accordingly, S^{q, cS) exhibits the spectral weight transfer, as 
seen in Fig. 1(b). 

We show S^{q) at /i = -1 and Ji = -2.5 in Figs. 4(b) and 
4(c), respectively. We find a sharp peak structure in the VC 
phase at small m, while the peak structure becomes broad in 
the SN/SDW phase at large m. This is because the transverse 
spin correlation exhibits an algebraic decay in the VC phase, 
and it is short-ranged in the SN/SDW phase. As m increases, 
the borad peak is leveled off, while the baseline goes up, and 
eventually we find a completely fiat profile, i.e., S^{q) = 1/4, 
at the saturation. On the other hand, the peak position depends 
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Eig. 5. (Color online) The extrapolated gap A(m) for several values of Ji 
at 72 = 1 as a function of m. Here we plot data in the SN/SDW phase. Note 
that the zero-field gap is tiny,^^^ and the VC phase is supposed to be gapless. 


on Ji and m. Indeed, we find a peak near q = nj2 regardless 
of m for = -1 [Fig. 4(b)], while a peak is at ^ = 25;:/64 
at m = 0 and it moves to ^ = 7r/2 near the saturation for 
/i = -2.5 [Fig. 4(c)]. 

To clarify the m dependence of the dominant q component 
of the transverse spin correlation, we determine a pitch angle 
Q from the maximum point of S~{q), as shown in Fig. 4(d). 
Note that S~{q) and S^{q) have a peak at the same position. 
In the VC phase at small m, Q agrees with the lowest-energy 
peak position q^ of S~{q, cS), as denoted by solid symbols in 
Figs. 3(b) and 4(d). However, in the SN/SDW phase at large 
m, Q behaves quite differently from q^ as m approaches the 
saturation. We find that Q changes toward 7r/2 regardless of 
/i, rather than the classical pitch angle cos“^(-/i/4/2) as q^, 
indicating that the dominant q component is not equivalent to 
the momentum position of the opening gap. This discrepancy 
is naturally understood in terms of the short-range nature of 
the transverse spin correlation. It is inadequate to take the only 
leading asymptotic term that represents the exponential decay 
of the transverse spin correlation function in order to describe 
the excitation dynamics correctly. In fact, we have seen that 
S^{q) turns into the fiat profile at the saturation, meaning that 
not only the dominant q component but also other components 
equally contribute to the spectral weight. 

Finally, we present the dependence of the energy gap of the 
transverse spin excitation in Fig. 5. We extrapolate finite-size 
data to the thermodynamic limit V ^ oo by assuming a linear 
relation A(N,M) = A(m) -f a/N. For m > 0.1, we find that 
A(m) increases as |/i| becomes large for Ji = -0.5, -1, and 
-1.5, indicating that the ferromagnetic exchange interaction 
stabilizes the two-magnon bound state. A(m) turns to decrease 
with further increasing |/i|, shown for Ji = -2. For small m, 
A(m) is rapidly reduced as m decreases down to the boundary 
between the SN/SDW and VC phases. 

In summary, we have studied the spin excitation dynamics 
of the frustrated ferromagnetic chain in the magnetic field by 
numerical methods. In the field-induced SN regime, the spin 
excitation spectra exhibit highly anisotropic behavior between 
gapless longitudinal and gapped transverse compoments. The 
field dependence of the gapless point of S^(q, oj) is consistent 
with the dominant longitudinal spin correlation. In contrast, 
the gap position of S^(q, oj) shows a unique field dependence 
that contradicts with the dominant transverse spin correlation. 
We hope that these features could be examined by inelastic 


neutron scattering experiments when searching for signatures 
of the SN/SDW phase. On the other hand, we should observe 
gapless excitations in the quadrupole channel. It would be an 
interesting future problem to study the quadrupole excitation 
dynamics and its possible relevance to observables. 
Acknowledgment The author thanks T. Masuda, S. Maekawa, M. Mori, 
and T. Sugimoto for discussions and comments. Part of computations were 
carried out on the supercomputers at the Japan Atomic Energy Agency and 
the Institute for Solid State Physics, the University of Tokyo. 


1) A. F. Andreev and I. A. Grishchuk, Sov. Phys. JETP 60, 267 (1984). 

2) For a recent review, see K. Penc and A. M. Lauchli, in Introduction to 
Frustrated Magnetism, ed. C. Lacroix, P. Mendels, and F. Mila (Springer- 
Verlag, Berlin, 2011), p. 331. 

3) A. V. Chubukov, Phys. Rev. B 44, 4693 (1991). 

4) L. Kecke, T. Momoi, and A. Furusaki, Phys. Rev. B 76, 060407(R) 
(2007). 

5) T. Vekua, A. Honecker, H.-J. Mikeska, and F. Heidrich-Meisner, Phys. 
Rev. B 76, 174420 (2007). 

6) T. Hikihara, L. Kecke, T. Momoi, and A. Furusaki, Phys. Rev. B 78, 
144404 (2008). 

7) J. Sudan, A. Ltischer, and A. M. Lauchli, Phys. Rev. B 80, 140402(R) 
(2009). 

8) M. E. Zhitomirsky and H. Tsunetsugu, Europhys. Lett. 92, 37001 (2010). 

9) M. Sato, T. Hikihara, and T. Momoi, Phys. Rev. Lett. 110, 077206 (2013). 

10) N. Shannon, T. Momoi, and P. Sindzingre, Phys. Rev. Lett. 96, 027213 
(2006). 

11) H. Tsunetsugu and M. Arikawa, J. Phys. Soc. Jpn. 75, 083701 (2006). 

12) A. Lauchli, F. Mila, and K. Penc, Phys. Rev. Lett. 97, 087205 (2006). 

13) T. Momoi, P. Sindzingre, and K. Kubo, Phys. Rev. Lett. 108, 057206 

( 2012 ). 

14) M. Enderle, B. Fak, H.-J. Mikeska, R. K. Kremer, A. Prokofiev, and W. 
Assmus, Phys. Rev. Lett. 104, 237207 (2010). 

15) T. Masuda, M. Hagihala, Y. Kondoh, K. Kaneko, and N. Metoki, J. Phys. 
Soc. Jpn. 80, 113705 (2011). 

16) L. E. Svistov, T. Fujita, H. Yamaguchi, S. Kimura, K. Omura, A. 
Prokofiev, A. I. Smirnov, Z. Honda, and M. Hagiwara, JETP Lett. 93, 
21 ( 2011 ). 

17) M. Mourigal, M. Enderle, B. Fak, R. K. Kremer, J. M. Law, A. Schnei- 
dewind, A. Hiess, and A. Prokofiev, Phys. Rev. Lett. 109, 027203 (2012). 

18) K. Nawa, M. Takigawa, M. Yoshida, and K. Yoshimura, J. Phys. Soc. Jpn. 
82, 094709 (2013). 

19) M. Hase, H. Kuroe, K. Ozawa, O. Suzuki, H. Kitazawa, G. Kido, and T. 
Sekine, Phys. Rev. B 70, 104426 (2004). 

20) A. A. Bush, V. N. Glazkov, M. Hagiwara, T. Kashiwagi, S. Kimura, K. 
Omura, L. A. Prozorova, L. E. Svistov, A. M. Vasiliev, and A. Zheludev, 
Phys. Rev. B 85, 054421 (2012). 

21) B. Willenberg, M. Schapers, K. C. Rule, S. Stillow, M. Reehuis, H. Ryll, 
B. Klemke, K. Kiefer, W. Schottenhamel, B. Buchner, B. Ouladdiaf, M. 
Uhlarz, R. Beyer, J. Wosnitza, and A. U. B. Wolter, Phys. Rev. Lett. 108, 
117202 (2012). 

22) S. E. Dutton, M. Kumar, M. Mourigal, Z. G. Soos, J.-J. Wen, C. L. Bro- 
holm, N. H. Andersen, Q. Huang, M. Zbiri, R. Toft-Petersen, and R. J. 
Cava, Phys. Rev. Lett. 108, 187206 (2012). 

23) K. Nawa, Y. Okamoto, A. Matsuo, K. Kindo, Y. Kitahara, S. Yoshida, S. 
Ikeda, S. Hara, T. Sakurai, S. Okubo, H. Ohta, and Z. Hiroi, J. Phys. Soc. 
Jpn. 83, 103702 (2014). 

24) M. Sato, T. Momoi, and A. Furusaki, Phys. Rev. B 79, 060406(R) (2009). 

25) M. Sato, T. Hikihara, and T. Momoi, Phys. Rev. B 83, 064405 (2011). 

26) O. A. Starykh and L. Balents, Phys. Rev. B 89, 104407 (2014). 

27) H. Onishi, J. Phys.: Conf. Ser. 592, 012109 (2015). 

28) A. Smerald, H. T. Ueda, and N. Shannon, Phys. Rev. B 91, 174402 
(2015). 

29) S. R. White, Phys. Rev. Lett. 69, 2863 (1992). 

30) E. Jeckelmann, Phys. Rev. B 66, 045114 (2002). 

31) J. des Cloizeaux and J. J. Pearson, Phys. Rev. 128, 2131 (1962). 

32) J. Ren and J. Sirker, Phys. Rev. B 85, 140410(R) (2012). 

33) C. Itoi and S. Qin, Phys. Rev. B 63, 224423 (2001). 

34) A. V611 and S. Wessel, Phys. Rev. B 91, 165128 (2015). 





